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Abstract 

This paper deals with a general method for the reduction of quantum systems 
with symmetry. For a Riemannian manifold M admitting a compact Lie 
group G as an isometry group, the quotient space Q = M/G is not a smooth 
manifold in general but stratified into a collection of smooth manifolds of 
various dimensions. If the action of the compact group G is free, M is made 
into a principal fiber bundle with structure group G. In this case, reduced 
quantum systems are set up as quantum systems on the associated vector 
bundles over Q = M/G. This idea of reduction fails, if the action of G on M 
is not free. However, the Peter- Weyl theorem works well for reducing quantum 
systems on M. When applied to the space of wave functions on M, the Peter- 
Weyl theorem provides the decomposition of the space of wave functions into 
spaces of equivariant functions on M, which are interpreted as Hilbert spaces 
for reduced quantum systems on Q. The concept of connection on a principal 
fiber bundle is generalized to be defined well on the stratified manifold M. 
Then the reduced Laplacian is well defined as a self-adjoint operator with 
the boundary conditions on singular sets of lower dimensions. Application to 
quantum molecular mechanics is also discussed in detail. In fact, the reduction 
of quantum systems studied in this paper stems from molecular mechanics. If 
one wishes to consider the molecule which is allowed to lie in a line when it is 
in motion, the reduction method presented in this paper works well. 
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I. Introduction 



Symmetry has always played an important role in mechanics. When a Hamiltonian system admits a 
symmetry group, the system reduces to a Hamiltonian system of less degrees of freedom. Reduction 
of Hamiltonian systems with symmetry is established by Marsden and Weinstein|l| and is explained 
in a textbook^]. The reduction method has found a wide variety of applications. In particular, 
the notion of moment map, which is a key word in the reduction method, has been introduced in 
differential geometry. 

The reduction method is applicable to molecular mechanics, both classical and quantum. In fact, 
the translational invariance allows one to separate the relative motion of atoms from the motion of 
the center of mass. Contrary to this, the rotational invariance of molecules is not able to separate the 
vibrational motion from the rotational motion of the whole molecular system, which was shown by 
Guichardetj|]. The nonseparability of vibration and rotation of a molecule is an underlying principle 
that allows cats to fall on their legs when launched in the air. While they have zero angular momentum 
and is free from external torque during the fall, they can make a rotation after a vibrational motion. 
It was Guichardet H who gave a decisive answer to this seemingly strange fact. He defined a natural 
connection on the center-of-mass system and thereby applied the holonomy theorem to show that 
the existence of nonvanishing curvature results in the nonseparability of vibration and rotation. It is 
to be noted here that the center-of-mass system is made into a fiber bundle if collinear configurations 
of the molecule are gotten rid of. The connection is defined on this restricted center-of-mass system. 

In spite of the nonseparability of vibration and rotation, the reduction method is still applicable. 
One of the authors (T.I.) applied the reduction method of Marsden and Weinstein to the molecular 
system by the use of conservation of angular momentum and gave a reduced Hamiltonian system [Q], 
according to which the "internal" motion of molecules is coupled with rotation through a kind of 
gauge field. The expression of the reduced equation of motion in terms of local coordinates is given 
in ||. He also constructed the Schrodinger equation to describe the internal motion of molecules in 
both two and three dimensions, using differential geometry of fiber bundles and connections^, 0, ^. 
However, in the previous papers || [7], ||, collinear configurations of the molecule in which all the atoms 
are aligned in a line in R? were out of consideration for the reason that the collinear configurations 
form a singular point set which prevents the center-of-mass system from being made into a principal 
fiber bundle. The bundle picture is extended to be applicable to a system of rigid bodies. Classical 
and quantum mechanics for the system of jointed identical axially symmetric cylinders is treated in 



It is Kummer[10| who first discussed the reduction of the cotangent bundle of a principal fiber 
bundle with a connection. He may have been aware of a use of the connection in many-body systems. 
However, he did not refer to the center-of-mass system as a principal fiber bundle in general, but 
referred to a planar three-body system as an S^-bundle. He claimed also that the reduction method 



had been applied to the planar three-body system by Satzer[ll]. Guichardet is the first who showed 
that the connection defined on the restricted center-of-mass system as a principal fiber bundle played 
an important role in the study of many body systems; he showed that the holonomy group of the 
connection coincides with the structure group by the use of the Ambrose-Singer holonomy theorem, 
along with the conclusion that any rotation can be realized as a holonomy associated with a closed 
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loop in the base manifold (or a shape space). 

This paper has an aim to improve the previous theory so that it may be applicable on the whole 
center-of-mass system which includes the singular point set stated above. Since the center-of-mass 
system is endowed with a natural Riemannian metric and admits the action of 50(3) as an isometry 
group, a general setting to start with is that a configuration space M is a Riemannian manifold on 
which a compact Lie group G acts by isometry. For a quantum system on M with symmetry Lie 
group G, the reduced quantum system by symmetry is to be defined on the quotient space Q = Mj G. 
However, a problem arises since the Q is not a smooth manifold in general, which may include singular 
points. A part of Q, denoted by Q^, is a smooth manifold, which is called the internal or shape space 
endowed with local coordinates describing the internal degrees of freedom of the molecule. One can 
set up a reduced quantum system on the smooth manifold Q^H 0- However, if one considers the 
whole Q, then a question arises as to how wave functions of internal coordinates should behave on 
the singular point set. This article will provide a general formulation to describe quantum mechanics 
of a reduced system with singular points taken into account. 

Let L^{M) be the Hilbert space of square integrable functions on the Riemannian manifold M, 
which is to be looked upon as the space of wave functions on the center-of-mass system. The group 
action G x M — > M; (g,x) \— * gx induces a unitary representation of G in L,2{M) through the 
action G x Lz(M) — > Lz(M); (g,f(x)) t— > f(g x). This representation will be decomposed into 
irreducible ones to give rise to representation subspaces of L/2(M) accordingly. To get an idea of the 
decomposition of L,2(M), the Peter- Weyl theorem on unitary representations of compact groups is 
of great help, since the theorem is understood to provide the decomposition of the Hilbert spaces of 
functions on groups. The decomposition of I^^M) is then performed in a similar manner, to define 
the spaces of equivariant functions on M according to respective irreducible unitary representations 
of G. These spaces will give rise to wave functions reduced by the use of the angular momentum 
conservation, since choosing an irreducible unitary representation of the isometry group amounts to 
keeping an angular momentum eigenvalue fixed. The equivariant functions will be shown to satisfy a 
good boundary condition, in a natural manner, on the set of singular points. The general formalism 
for reducing wave functions is applicable to molecular mechanics. In particular, a triatomic molecule 
will be studied in detail. 

The idea of using the Peter- Weyl theorem for quantization is old; it dates back to Casimir's 
work on quantum mechanics of rigid body in 1931. The idea of using the invariant functions to 
set up a reduced quantum system by symmetry is also old; for example, the quantum theories of 
electromagnetism[12] or gravitation [13, 14] are described by gauge invariant and reparametrization 
invariant states, respectively. The idea of using equivariant functions is not new, either; it has already 
been introduced by T.I.||, |7j in molecular quantum mechanics and by Landsman and Linden [|15[] in 
quantization on homogeneous spaces, respectively. A point to make in the present paper is to extend 
the concept of equivariant functions on a principal fiber bundle to any manifold on which a compact 
group dose not necessarily act freely. 

There are other methods of reduction by symmetry. A more algebraic formulation was set up by 
Landsman [l(| and Wren[[l7]], who used a representation theory of C* algebras due to Rieffelf^]. The 
present formulation is rather geometric and comprehensible in applications. On the other hand, the 
path integral formulation is compatible with a method of reduction by symmetry. It has been shown 
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Jig , |i"9| , [20[ | that path integral over a Lie group G is reduced to path integral over a homogeneous space 
Q = G/H when the system possesses symmetry given by a subgroup H. However, this method is 
only applicable to homogeneous spaces, which are free from singularity, while the method proposed 
in the present paper will be applicable even when the quotient space Q = M /G is not a smooth 
manifold, as will be shown later. 

On the other hand, Sniatycki and Weinstein[21], and other people[22] have studied reduction and 
quantization of symplectic systems with symmetry. They consider reduction of symplectic manifold 
when the momentum map takes a singular value and therefore the level set does not form a smooth 
manifold. In our context, singularity refers to points of a configuration space which admit larger 
isotropy groups than those at generic points. The subject to be considered here is different from 



what has been considered in 21], and a relation between these subjects will not be examined. 



The plan of the present paper is as follows: In Sec. II, the Peter- Weyl theorem on unitary 
representations of compact groups is reviewed briefly. This theorem is extended to be applied to 
a unitary representation of a compact Lie group G in a Hilbert space TC in order to decompose TC 
into a series of invariant subspaces. The decomposition procedure is then applied to the unitary 
representation of G in the Hilbert space I/2(M). As far as the decomposition is concerned, M is 
assumed only to be equipped with a G-invariant measure, and does not need to be a Riemannian 
manifold. The above stated decomposition of L2(M) proves to be a decomposition into a series 
of spaces of "equivariant" functions. Section III contains a study of geometric structure of the G- 
manifold M, where the assumption is not yet made that M is a Riemannian manifold. It will be 
shown that the orbit space Q = M/G becomes not a smooth manifold but a collection of smooth 
manifolds of various dimensions, which structure is called stratification. With this stratification 
taken into account, a connection on M , equivariant functions on M, and covariant derivatives will be 
defined. In Sec. IV, M is assumed to be a Riemannian manifold, and G to act on M by isometry. The 
Laplacian acting on smooth functions on M will be reduced to be defined on each space of equivariant 
functions according to the decomposition of L2 (M) , which will provide a reduced quantum system on 
Q. The reduced systems are well defined. In fact, the equivalence condition will provide a boundary 
condition on equivariant functions at singular points. As a simple example, the Laplacian on R 2 
will be studied. In Sec. V, the general formalism developed in the preceding sections is applied 
to molecular mechanics. A general setting for iV-atomic molecules will be established, and then 
triatomic molecules will be studied in detail. Sec. VI contains conclusion and discussion. 



II. Reduction of quantum systems with symmetry 
A. The Peter- Weyl theorem 

We start with a brief review of the Peter- Weyl theorem, which will provide a key idea to the reduction 
method for quantum dynamical systems with symmetry. 

Let G be a compact Lie group. Then there is the unique Haar measure [Xq on G which is 
normalized to satisfy f G dfj>G = 1- Let -^(G) denote the space of all the square integrable functions 
on G, L2(G) := {/ : G — > C | J G \ f{g)\ 2 d^G{g) < which is equipped with the inner product 

(/,/'):= / W)f'(gW G {g) (2.1) 

JG 
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for / and f of L 2 (G). 

Let p x : G — > U(TC X ) denote an irreducible unitary representation of G on a Hilbert space 
TC X of finite dimension cf x := dimTC x , where \ is to label all of inequivalent irreducible unitary 
representations of G, and U(TC X ) denotes the space of unitary operators on Tt x . By pfj(g) we denote 
the matrix elements of p x (g) with respect to some orthonormal basis of TC X , where indices range over 
i,j = l,2,-..,cP<. 



The Peter- Weyl theorem then states that all the functions {\/# P?j(g)}x,i,j f° rm a complete 
orthonormal set (CONS) in L 2 (G). Namely, one has the orthonormality relations 



d x / pfj&pUgWaig) = 6 XX 5 lk 5 j£ (2.2) 

J G 



along with the completeness condition that if 



pU9)f(9)dp G (g) = (2.3) 



G 11 



for all XihJi then / = 0. Hence, any function / E L 2 (G) can be expanded in a Fourier series in 
{y/dX pfj} x ,i,j, so that one has 

f(h) = Y. dX p¥ h )j WdfiaWaig) 
= E dX [ Pl(hg- l )f{g)dp G {g) 

= J2 dX [ pl^fig-'hWdg). (2.4) 



The expansion formula Q2.4| ) is put formally in a compact form, 

J2d x pU 9 ) = S(g), (2.5) 

xA 

where 5(g) is Dirac's delta function on G with respect to the measure pc- The Peter- Weyl theorem 
also implies that 

L 2 (G)^®((H X T®H X ), (2.6) 

x 

where it is to be noted that (TC X )* <S> T~C X is isomorphic with the direct sum of d x copies of TC X . 
B. Method of reduction 

We now apply the Peter- Weyl theorem to a quantum dynamical system with symmetry to obtain 
a series of reduced systems. A quantum dynamical system is defined to be a pair (TC, H) of a 
Hilbert space TC and a Hamiltonian H, where if is a self-adjoint operator on TC. Suppose that the 
system (TC, H) admits a compact Lie group G as a symmetry group, namely, each element g £ G is 
represented as a unitary operator U(g) which acts on TC, commuting with H. 
In view of ( |2.4| ), we define, for each label (x,i), & n operator P x on TC to be 

i?:=d* / pl(g)U(g)dp G (g), (2.7) 

JG 
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which apparently commutes with H. From the G-invariance of the measure pa an d from the or- 
thonormality relations ( |2.2| ), it follows that 



(if )t = px f pxpx' = §X x' S .. P x (2.8) 
Further, the completeness condition ( |2.5D implies that 

£^ x = id w , (2-9) 



where id^ is the identity operator on Ti. Equations (2.8) and (|2.S|) show that the set {P*} x ,i forms 



a family of orthogonal projection operators on Ti, bringing about the orthogonal decomposition of 

n, 

H = ®Imif, (2.10) 

x,i 

which is an analog to ( |2.6| ). Moreover, each subspace Im P x is invariant under the action of the Hamil- 
tonian H. Thus the dynamical system (Ti,H) is broken up into a family of subsystems (Im P*,H) 
labeled by (x,i)- We call each system (Im P*,H) a reduced quantum dynamical system. 

In the language of physics, the pair (x, i) labels conserved quantities associated with the symmetry 
group G, and thereby define a closed dynamical system that consists of the states with the assigned 
conserved quantities. For example, if the original system has G = SU(2) symmetry, the angular 
momentum is conserved. The states labeled by (j, m) have the total angular momentum and the 
component of the angular momentum fixed at J 2 = j(j + 1) and J3 = m, respectively, and are 
described as vector-valued functions with 2j + 1 components. 

C. Characterization of the reduced system 

To gain a deeper insight into the decomposition ( 2. ICQ , we introduce an operator on Ti by 



V*:=d* / pf J (g)U(g)d^ G {g) (2.11) 

J G 



for each label \ and indices i,j = l,---,d x . In particular, one has = P*. A straightforward 
calculation shows that 

(V0 = Vl V*V* = 8*6 jk VZ. (2.12) 



ijJ v ij v kl u ]k v i(- 

As an immediate consequence, one obtains V%P£ = if k 7^ j, so that the domain of V^f reduces 
naturally to ImPf on account of (pTTcj ). Since (VffiV* = VfiV* = V* = Pf, it holds that KerV$ = 
KexPf. Similarly, from Vfi(V0 = V$V% = V* = if , it follows that Ker {V^ = Ker if , and 
thereby that ImV*- = (Ker (V^) 1 = (Ker if ^ = Im(Pf)^ = ImPj X . Therefore, V*- becomes a 
unitary transformation 

Vj* : ImPf -► half. (2.13) 
As a collection of V?-, we define another operator V x : Ti x (g> Ti — ► Ti x <S> Ti by 

V x := d x [ p x (g) ® U{g)dn G (g). (2.14) 
JG 



Then Eq.fl2.r2Q implies that 



(yx)t = yx v x V x = d x V x , (2.15) 
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which shows that V x /d x is a projection operator on 7i x ® 7i. Further, it is easy to verify that V x 
satisfies, for any h £ G, 



V x (p x {h)®U(h)) = V x , (2.16) 

(p x (h) U(h))V x = V x , (2.17) 

which implies that 

lmV x = (H x ® Hf := {tp e U x % H | (p x (h) ® U{h))^ = ip, Wh e G}. (2.18) 

We call (Tt x (8> W) G the subspace of invariant vectors of H x <8> 

Let {e*, • • • ,e^ x } be an orthonormal basis of 7i x , which defines an injection ef : C — > W x by 
A i — > Ae* for each i = 1, •••,(#*. Its adjoint operator (e*)T : TC X —>■ C is defined by a surjection 
w i — ^ (e x ,v), v G The domain of respective maps extends to the tensor product space with 7i to 
give rise to ef : TL — ► TL X TL and (efy ■ TL X (8> TL — » With these notations, V^f is put in the form 
^ij = ( e f)^ ^ x e j • We then introduce an operator by 

Sf := —3= V x e x :TL ^TL X ®TL. (2.19) 

The adjoint operator is expressed as 

(S*)t := _L (e*)t yx : ft* ® ft _> ft. (2.20) 

VdX 



Then the second relation of ( |2.15| ) yields 

(SffSf = (ej)t e * = Pf. (2.21) 
On the other hand, from ( |2,12|) , we observe that V x ej (ej)^ V x = V x , and thereby obtain 

Sf(Sf^ = ^V x ej(e^V x = ^V x (2.22) 



Since Pf and V x /d x are projection operators on ft and on Tt x (g) ft, respectively, Eqs.( 2.21 ) and 
( 2.22| ) are put together to imply that the restricted map 

Sf : Im Pf AlmF = (H x TL) G (2.23) 

is a unitary transformation. Thus the subspace lm Pf is characterized as (TC X <8> TL) G . 

We turn to the Hamiltonian H acting on ft. It appears that H is extended to an operator 
id-ftx <8> H on TL X (g) ft. Since -ft commutes with U (g) for each g £ G, the extended operator id-^x ® H 
commutes with y x and with as well, so that one has 

SfH = (id„ x ® #)S* (2.24) 

It then turns out that the reduced system (lmP x , H) is identified with ((ft x (g)ft) , id^x ®H), which 
will be investigated in the following sections in detail. 

Before investigation into the reduced system, we wish to consider how the action of G is de- 
composed according to the decomposition of ft. In the below, the dual space (TL X )* to TL X is taken 
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as (ft x )* := Horn (ft x ; C), irrespective of the inner product in ft x . Then the map vd* (S*y = 
(e x )t V x :H X ®H -> H can be viewed as a map Tj : ft — > (ft x )* ft in the manner as follows: 
Using the identification End ft x = (ft x )* ® ft x , we may regard the map V x : ft x (£> ft — ► ft x ft as 
a map U x : ft -» (ft x )* ®H X ®H. Combining id (Wx) * ® (e£)t : (ft*)* ®H X ®H^ (ft x )* ft with 
V x , we can express T x as the map 

Tf = (id (wx) . ® (e*)t)V* : ft (ft x )* ® ft. 

Then it can be verified that 



id 



(2.25) 



(2.26) 



where (T x )t : (ftx) 
transformation 



• ft — > ft is the adjoint operator. From this, it follows that Tj yields a unitary 



Tf : Im P x ^(H x )* ®lmPf. 



(2.27) 



Then the right invariance of V x under the G action, expressed in ( |2.16| ), makes the following diagram 
commutative, 



ft = © Xii Imif — ^ (ft x )*®ImP x 



U(h) 



ft 



j 



|*p x ('i" 1 )«>id 



(2.28) 



(ft x )*<g)ImP/, 



where t p x (h 1 ) = p x (h) is the contragredient representation of G on (ft x )*. As for the action of the 
Hamiltonian H, we obtain the following commutative diagram accordingly, 



ft 

H 

ft 



(ft x )*®ImP> 



(2.29) 



X)< Imif 



The commutative diagrams ( |2.28| ) and ( 2.29 ) show that the decomposition of the representation 
(G, ft, £/) is compatible with the spectral resolution of H. 

D. Equivariant functions 

The general reduction method introduced in the previous section applies to a quantum system on a 
configuration space M which admits the action of a compact Lie group G. 

Suppose that a compact Lie group G acts on a differentiable manifold M by diffeomorphisms, 
namely, we are given a C°° map G x M — > M satisfying ex = x and (gh)x = g(hx) for any x € M and 
any g,h £ G with the identity element e £ G. Then M is called a G-manifold. In addition, suppose 
that M is equipped with a measure /ijn which is invariant under the action of G. Let ft = L,2(M) 
be the space of square integrable functions on M with respect to the measure fiM- The group G is 
represented on ft = L,2(M) by unitary operators £7(/i), h £ G, through 



(U(h)f)(x) :=/(rt), /€L 2 (M). 
8 



(2.30) 



Then the reduction method applies to Ti = L 2 (M) to yield the Hilbert subspace ImP x . The unitary 
transformation Sj given in ( gjp allows us to identify ImP* (c L 2 (M)) with (H x L 2 {M)) G (c 



TL X ®L 2 (M)). The space TL X ®L 2 (M) can be identified with the Hilbert space, L 2 (M;TC X ), of square 
integrable 7^ x -valued functions on M; 

L 2 (M;H X ):={^:M^H X \ [ U{x)\\ 2 d(i M (x) < oc}, (2.31) 

J M 

which is equipped with the inner product 

<0,V> := / (4>(x),^(x))dfi M (x), 0,4? e L 2 (M;H X ), (2.32) 



where (cj)(x),ip(x)) denotes the inner product in 7i x . Then, the condition given in ( [2.18 ) along with 



fl2~30D implies that ip G (H x ® L 2 (M)) G , when viewed as an "H x -valued function, satisfies 

^(hx) = p x (h)^(x), /iGG, (2.33) 

which shows that ip is equivariant under the G-action. We conclude therefore that the reduced 
Hilbert space ImP* is identified with the space of square integrable equivariant functions, which we 
denote by L 2 (M;H X ) G ; 

(W^L 2 (M)) G 9* L 2 {M;H X ) G 

(2.34) 

= & e L 2 (M;H X )\ iP(gx) = p x (g)iP(x), <? G G}. 



Here we have to note that according to the decomposition ( |2.10| ) along with ( |2.23j ), ( |2 . 2 7| ) , and ( [2.34 ). 
the Hilbert space L 2 (M) is decomposed into 

L 2 {M)^Q){(H X )* ®L 2 {M;H x f). (2.35) 

x 

So far we have characterized the reduced Hilbert space on the G-manifold M. We now have to 
specify the Hamiltonian on L 2 (M) and to reduce it. We will take a Hamiltonian as defined to be 
the sum of the Laplacian on M and a potential energy function on M. In the succeeding sections, 
we will study the geometric structure of M in order to analyze the Laplacian on M according to the 
decomposition ( |2.35| ). 

III. Geometric setting on G-manifolds 

A. Stratification of G-manifolds 



According to Davis[23], G-manifolds can be viewed as collections of fiber bundles. We here make a 



brief review of his idea in a suitable form for our application. For more rigorous definitions, see the 



literature [24]. 



Let M be a G-manifold. For a point x G M, we denote the isotropy subgroup at x and the G-orbit 
of x by G x := {g 6 G | gx = x] and by O x := {gx \ g G G}, respectively. Then one has O x = G/G x . 

Take G x and G y for two points x, y G M. If G x is conjugate to G y by an inner automorphism 
A g : G — ► G; h i— > ghg~ x with some g G G; G y = g • G x • g~ 1 , then the orbits O x and O y are 
diffeomorphic to each other by the correspondence induced by A g . Of course, if y is in the orbit 



O x , namely, if there exists some g G G such that y = gx, then G y is conjugate to G x by A g . The 
conjugacy class of G x is called the orbit type of O x and denoted by [G x ] or r. We say that the point 
x itself also has the orbit type r, if O x has the orbit type r. Let T(M) denote the set of all the orbit 
types appearing in M; T(M) := {[Gy | x G M}. For each r G T(M), we denote a representative of 
the conjugacy class by G T G r. 

One can introduce a partial order on T(M) as follows: We say that t\ is lower than T2 (tx,T2 G 
T(M)), ri < T2, if there are representatives G T1 G t\ and G T2 G T2 such that G T1 D G T2 . An orbit 
O x is called maximal if its orbit type r is maximal with respect to this order. Orbits which are 
not maximal are called singular. We say also that a point x is maximal or singular according to 
whether the orbit type of O x is maximal or singular. For each orbit type r G T(M), we denote by 
M T := {x G M I G x G r} the set of points with the same orbit type r, which becomes a smooth 
manifold. Thus M is stratified into a collection of the smooth submanifolds M T , M = TJ T M T , which 
is partially ordered by the system of orbit types (T(M), <). Each manifold M T is called a stratum. 

The set of orbits Q := M/G becomes a topological space with respect to the quotient topology 
which is defined by demanding that the canonical projection map tt : M — > M/G is continuous. The 
space Q is called an orbit space. A point q G Q is also classified by the orbit type of ■K~ 1 (q). The 
point q is said to be maximal or singular according to whether the orbit is maximal or singular. The 
Q inherits differentiable structure from M, if and only if all the orbits have the same orbit type. 
If otherwise, Q is stratified into a collection of smooth manifolds of various dimensions. Setting 
Q T := ir(M T ), one has Q = TJ r Q T . The restriction of tt : M — > Q onto each stratum n T : M T — > Q T 
defines a fiber bundle with fiber G/G T . We then call the pentaplet (M,G,Q,tt,T(M)) a stratified 
bundle. 

If every point of M has the same orbit type, tt : M — > Q is nothing but a usual fiber bundle. 
In particular, if the group G acts on M freely, namely, G x = {e} for all x G M, then tt : M — > Q 
becomes a principal fiber bundle with structure group G. In this sense, we may regard the G-manifold 
as a generalization of fiber bundles although the base space Q = M/G is not a smooth manifold. 
Mechanics of molecules provides an example of a stratified bundle that is not a fiber bundle, as shown 
in the later section. 



According to the principal orbit theorem] 25, 23], if the orbit space Q = M/G is connected, there 
is the maximum orbit type in T{M) with respect to the order <. Although the maximum orbit 
type is also named the principal orbit type in [25, 23 1, we will call it maximum in this paper. We 
assume that Q is connected and that the maximum orbit exists. We denote the maximum orbit 
type by /i. Moreover, the principal orbit theorem|25, ^] states that the maximum stratum is 
open and dense in M. Thus the set M s ; ng := M — M^ of all the singular points coincides with the 
boundary dM^. The image of M M and dM^ by the projection tt are denoted by := ir(M^) and 
dQn ■= 7r(<9M M ), respectively. We put dimM = m and dimQ^ = n = m — p. Then the dimension of 
the maximum orbit is dimG/G^ = p. 

Though the orbit space Q is not a manifold, one can speak of differentiability of functions on Q. 
A function of (p : Q — > R is called of C r class when (p o tt : M — > R is a differentiable function of C r 
class. Clearly, cp o n is a G-invariant function on M. Conversely, any G-invariant function / on M is 
identified with a function on Q. We denote the space of G-invariant functions on M of C r class by 
C r (M) G = {f:M -^R\f(gx) = f(x), V 5 G G, Vrc G M}. 
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A tangent vector X to M at x is usually defined as a differential operator acting on C°°(M); 
X : C°°{M) — > R. A tangent vector on Q is defined as follows: We define an equivalence relation 
~ in the tangent vector space T X M by stating that X ~ Y if Xf = Yf for all / G C°°(M) G . The 
equivalence class of X is denoted by n*(X), which defines a linear map tt*(X) : C°°(M) G — > i?. The 
set T g Q := {7r*(X) | X G T X M, 7r(x) = (?} becomes a vector space through the structure of the vector 
space T X M and is called a tangent vector space at q G Q. 

B. Stratified connection 

Let (M, G, Q, ir, T(M)) be the stratified bundle defined above. If G acts on M freely, the stratified 
bundle becomes a principal fiber bundle. Although the concept of connection is usually defined on 
principal fiber bundles, we would here like to define extended connections on stratified bundles. 

Let us call a subspace V x := T x O x of T X M a vertical subspace at x G M. The action of g G G : 
M — ► M; x \—* gx induces an action g* : T X M — > T gx M by differentiation. The complement H x 
of V x is called a horizontal subspace. If the direct sum decomposition T X M = V x © H x is smooth 
and if the family {H x } x ^m satisfies the invariance; H gx = g*H x , the decomposition T X M = V x © H x 
defines a connection. However, we should note that dim 14 = dimC4 and dimi^ = dimM — dimO x . 
jump suddenly when the point x passes singular points. Thus the smoothness of the decomposition 
T X M = V x ®H X must be required on each stratum M T , so that one understands that for any smooth 
vector field X which is decomposed into X (x) = X (x) + X H (x) according to T X M = V x © H x , the 
components X v and X H are also smooth on each M T . Then the assignment x i— > H x for each x G M 
is called a stratified connection over the stratified bundle ir : M — > Q. 

The decomposition T^M = 14® induces a decomposition of the dual space T*M = V* ® il* 
with I/* := {0 € T X *M | = 0, W G F K } and H* = {ip G T X *M | V(n) = 0, Vn G 14}, 

Let g and g x denote the Lie algebras of G and G x , respectively. The relation g gx = Ad^g^, is an 
immediate consequence of G gx = A g G x . The group action G x M — » M; (<7,x) ^ gx gives rise to 
vector fields g x M — > TM; (^, x) i— > £m( x ) as infinitesimal transformations. Fixing a point x G M, 
one obtains a linear map 9 X : g — > T X M; 9 X (£) = £m( x )- It then follows that Ker^ = g x and 
lm9 x = 14, and hence that 9 X : g ^ T X M induces an isomorphism 9 X : g/g x — * V x . 

The connection defined above is described in term of differential forms. A connection form 
u is defined as the composition of the projection T X M = V x © H x — * V x and the inverse map 
(9 x )- 1 :V x ^g/g x : 

u x :T X M = 14 © H x -> 14-^ g/g, (3.1) 

at each point x G M. The form to is thus a one-form which takes values in quotient spaces of the Lie 
algebra g. It has the following properties: 

u(9 x (0)=Z (modgj, £Gg, (3.2) 
(g*u) x = Ad g uj x , g G G, (3.3) 

where g* is the pull-back associated with the map g* : T X M — ► T gx M and Ad 9 is to be understood 
as a map Ad 9 : g/g 2 — > g/g gx . To verify (3.3), we need the formula that 9 gx (Ad g ^) = g*(9 x (£)) for 
5 £ g and that (^x) -1 9* = Ad 9 (0 a; )~ 1 . The properties ( |3.2[) and (|3.3| ) are generalization of the 
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well-known denning properties of usual connection forms. It is also noted that the composition map 

T x o uj x : T X M = V X ®H X ^V X (3.4) 
is a projection, and that Eq.(3.2) is equivalently written as u x o 6 X = idg/g . 

C. Equivariant forms 

Let tp G C°°(M,7i x ) G be a smooth equivariant function, that is, one has t^(gx) = p x (g)^(x) for 
g G G (see (fE33p). If a point x G M carries a nontrivial isotropy group G^, then the value of ip(x) 
becomes invariant under the action of G x ; ip(x) = p x (g)ip(x) for g G G x . For a subgroup Gi C G, 
we here define (7i x ) Gl to be a maximum subspace of invariant vectors under the action of Gi, that 
is, (H x ) Gl := {v G H x | p x (c/)f = v,\/g £ G x }. With this notation, we then have ^(z) G {H x ) Gx . 

Properties of subspaces of invariant vectors are worth remarking. One has obviously (TC X )^ = 
7i x . Further, one obtains (7i x ) G = {0}, if and only if there is no nontrivial invariant vector. To a se- 
quence of subgroups {e} C G\ C G^C G), there corresponds a sequence of subspaces 7i x D (W*) Gl D 
(H X ) G2 (D {0}). For conjugate subgroups Gi and A,Gi, it holds that (H x ) AgGl = p x (g)(H x ) Gl . In 
particular, on an orbit O x , one has (H X ) G ^ = p x {g){H x ) G * on account of G gx = A g G x . 

Like equivariant functions, we can define equivariant differential forms in a similar manner. An 
equivariant k-form is defined as an 7i x -valued differential /c-form a : f\ k TM — > W x satisfying g*a = 
p x (g)a for any g £ G. When an equivariant fc-form a further satisfies i(^A/)a = for any £ G g, 
it is called an equivariant horizontal k-form. While the set of all the differential /c-form on M is 
denoted by f2 k (M), the set of all the W x -valued differential /c-forms and the set of all the equivariant 
horizontal A:-forms are denoted by H x ® Q k {M) ^ f2 k (M;Ti x ) and Q k H {M] H X ) G , respectively. 

At each point x G M carrying a nontrivial isotropy group G x , the equivariant horizontal A;-form 
a takes a restricted range as well as equivalent functions. It turns out that a x {X\, • • • , X^) G (7i x ) Gx 
for any X\, • • • , X^ G T X M when G x is connected. The proof is given below: Let Lx denote the Lie 
derivation by a vector field X on M, and p x denote a representation of the Lie algebra g on TC X 
induced by differentiation of the representation p x of G. The defining property of the equivariant 
form g*a = p x (g)a for g G G is differentiated to give L^ M a = p x {C)a for £ G g. This equation and 
the horizontal condition i{ihi)a = for £ G g are put together with the Cartan formula L^ M a = 
i(£,M)da + di(^M)ct to provide 

P*(0 a = i(^M)da, £ G g. (3.5) 

If £ G gj,, one has ^m{x) = and hence p*(^)a x = from (|3,5|), which implies that p x (g)a x = a x for 
g G G x if G x is connected. Thus we conclude that a x G (7i x ) x . 

D. Associated vector bundles 

So far we have discussed equivariant functions and forms on the stratified bundle. We now wish 
to define vector bundles associated with the stratified bundle, like vector bundles associated with 
principal fiber bundles. 

Let us define an equivalence relation ~ in M x TC X by the relation (x,v) ~ (gx, p x (g)v) with 
g G G. Let [x,v] denote the equivalence class with a representative (x,v). Then the vector bundle 
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associated with the stratified bundle (M,G,Q,tt,T(M)) by a representation (7i x ,p x ) is defined to 
be the quotient space 

M x pX H* := ( I] (to x (n x ) G *)) I ~ • (3.6) 

The projection map M x W x — > M naturally induces a projection map 7r p x : M x pX TL X — > Q; [x,v] i— ► 
7r(x). Further, each point x G -M defines an isomorphism x : (7i x ) Gx — > ir p x (ir(x)) = O x x pX H x by 
v i — y x(v) := [x,v]. Note that, for each stratum Q T , TTpxiQr) is a vector bundle over Q T with fiber 
(TC x ) Gt , so that one has M x pX W x = U T 7r^x (Qt)- I n this sense, we may call M x pX 7i x a stratified 
vector bundle. However, we will refer to it as the associated vector bundle for simplicity. Moreover, 
each fiber ft^x(q) at q G Q inherits an inner product from 7i x ; for 77,77' G TT~x(q), the inner product 
(77, 77') is defined to be 

(77,77') :=(x- 1 (77),i- 1 (77')), (3.7) 

where the RHS is the inner product defined on TC X . It is easy to verify that the RHS is independent 
of the choice of x G 7r _1 (g). 

A section of the associated vector bundle is a map a : Q — > M x pX TL X satisfying ir p x o a = idg. 
An equivariant function ip defines a section 1$ of the associated vector bundle through ip b (ir(x)) = 
[x, i/j(x)\ = (xoip)(x). The ip (q) is well-defined. In fact, ip b (Tr(gx)) = [gx,ij)(gx)] = [gx, p x (g)^p(x)] = 
[x,i/)(x)\ = ^(tt(x)) for g G G. Conversely, a section a defines an equivariant function cr* through 
<t*(x) = (x _1 o a o 7r)(x). Thus we can identify sections of the associated vector bundle with 
equivariant functions. When the corresponding equivariant function cr* is differentiable as an TL X - 
valued function on M, the section a is called differentiable, while Q is not a manifold in general. 
The set of all of the differentiable sections of M x pX TC X is denoted by r(M x pX TC X ). Thus the space 
of sections r(M x pX TC X ) is in one-to-one correspondence with the space of equivariant functions 
C°°(M;n x ) G . 

To associate equivariant horizontal fc-forms with certain sections, we should set up another vector 
bundle. The action of g G G on the space of one- forms, (g- 1 )* : T*M -> T* X M, induces the action 
on the space of horizontal fc-forms, (g^ 1 )* ■ f\ k H* —* A k H* x . This action is extended to that on the 
space of W x -valued horizontal fc-forms A k H* (g) TC X by the linear map generated by 

(j) ® v G A k H* x ®H X ^ p%)(</> ® v) := {g' l )*<t) p x (g)v G A k H* x H x . (3.8) 

We then take a subspace of invariant vectors 

(A k H* x ® H x f* := {C G A k H* x ® H x \ p*(g)( = (, Vg G G x } (3.9) 

and define an associated vector bundle, like ([[!]), through 

M x~ (A k H* ® H x ) := ( [] ({x} x (A k H* x ® W*) G * J / ~, (3.10) 

where the equivalent relation ~ in M x (A k H* ®7i x ) is defined as (x, Q ~ (gx, p x (g)Q- A projection 
map 7r^ : M x ~ (A k H* ® H*) -> Q is defined naturally as [x, C] l— ► 7r(x). The space of smooth 
sections r(M x~ x (A k H* (£> T~i x )) is also in one-to-one correspondence with the space of equivariant 
horizontal forms £2jj(M; Tt x ) . 
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E. Covariant derivative 

Covariant derivatives of fc-forms on the stratified bundle can be denned like those on the principal 
fiber bundle. Covariant derivation is a linear map D : Q k H {M-H x ) G -> Q k H +1 (M;H x f defined 
through 

Da(X 1 ,---,X k+1 ) = da(X?,---,X k H +1 ), X t E T X M, i = 1, • • • , k + 1, (3.11) 

for a E f2fl(M;?i x ) G , where is the horizontal component of X% = Xj + Xf . 

The covariant derivative Da can be expressed by using the connection form u. To any X v E V x , 
there corresponds an element £ E g such that X y = uniquely modulo g^. By the definition of 

the connection form u>, we then have £ = uj{X v ) = lo(X) (modgj.). With the help of ( |3.5| ), one has 
i(X v )da = i(£,M)da = p x (£)a = p x (ui(X))a. The last equality holds well in spite of the ambiguity 
in the value of u(X) since p x (C)a x = for £ E g x . Putting together these equations results in 

Da(X u ---,X k+1 ) 

= da(X\ — X x , ■ ■ ■ ,X k+1 — X k+ i) 

fe+i 

= da(X lr --,X k+1 ) + Y / (-l) l P X (u(X i ))a(X u ---,X i ,---,X k+1 ), (3.12) 

8=1 

where Xi means the removing of Xi . Thus our result is expressed as 

Da = da — p x {oj) A a. (3.13) 
We should note that smoothness of the covariant derivative is ensured only within each stratum M T . 

IV. Reduction of the Laplacian 

A. Reduced Laplacians 

We have characterized the Hilbert space of the reduced quantum system and studied the geometric 
structure of the manifold which is underlying the Hilbert space. Now we turn our attention to the 
Hamiltonian acting on the Hilbert space Li2(M). As was anticipated, the Hamiltonian H we treat 
takes the form H = i^^-m + V with the Laplacian Am and the potential V. 

First we give a precise definition of the Laplacian. We have been working with the G-manifold M 
equipped with the G-invariant measure [Im- We assume that M has no boundary. In what follows, 
we make another assumption on M; M is assumed to be oriented and equipped with a Riemannian 
metric qm which is invariant under the action of G. The metric gu induces a volume form x>Mi which 
is also invariant under the action of G. We assume also that pm is the measure associated with the 
volume form vm- Let C^°(M) be the set of all the C°° functions on M with compact support. Then 
the Laplacian A M : C C °°(M) -> C™(M) is defined through 

/ \\df{x)\\ 2 9M v M = [ f(x)(A M f)(x)v M , (4.1) 

JM y JM 

where ||4f(^)||g M denotes the norm of T*M induced by the metric gu- We note that the domain 
of Am is extended in L,2(M) to make Am a self-adjoint operator. Of course, in order that this be 
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the case, M has to be assumed to be complete. Since both the metric gM and the volume vm axe 
G-invariant, the Laplacian Am is also G-invariant; namely, Am commutes with U{g) for any g £ G. 

Next we turn our attention to the potential energy V(x). It is a smooth function V : M — > R 
acting on / G L 2 (M) by multiplication as (Vf)(x) := V(x)f(x). We assume that V is also G- 
invariant; V{gx) = V{x) for any g £ G and x £ M. Thus the action of V also commutes with 

U(g). i 

Since each term of the Hamiltonian H = — Am + V commutes with U(g) for any g £ G, we can ap- 
ply the decomposition ( [2.10D and the commutativity ( [2,24 ) to both Am and V separately. Hence, Am 



and V act as (id W x©A M ) and (id H x®V) on the reduced Hilbert space {H X ®L 2 (M)) G = L 2 (M; H X ) G . 
The Laplacian (\d H x ® A M ) with the domain restricted to (H X ®L 2 (M)) G = L 2 (M;H X ) G is denoted 
by A x , and is called a reduced Laplacian. The reduced Laplacian A x is also a self- adjoint operator. 
Then for an equivariant function ip £ C£°(M; 7i x ) G , the defining equation of the Laplacian A x takes 
the form 

" \\d^)\\L r VM= [ mx),(^)(x))vM, (4.2) 



M 9M JM 



where in the LHS ||c?'0( a; )llg M denotes the norm of TL X (&T*M induced from the metric gM- 
B. Rotational and vibrational energy operators 

To make further study of A x , we make intensive use of the vertical-horizontal decomposition T X M = 
V x © H x introduced in Sec. III. We have not chosen a specific connection yet. Now we fix the 
connection by demanding the orthogonality V x _L H x with respect to the metric gM(x)- Then, the 
G-invariance of the metric ensures that g*H x = H gx , and hence a unique connection is determined. 
Since the set of maximum points is an open dense subset of M as noticed previously and since 
the connection is smooth when restricted to M„, the Laplacian on will be smoothly decomposed 
into two, vertical and horizontal components, by the use of the connection. 



According to the orthogonal decomposition T*M = H* © V*, the integrand of the LHS of (4.2) 
is also written as 

\\di>(x)\\ 2 gM = \\mx)) H \\ 2 gM + \\mx))v\\ 2 9M . (4.3) 

By the definition of the covariant derivative (|3.11] ) and the expression ( |3.13| ) in terms of the connection 
form, the above equation is put in the form 



\mx)\\ 2 9M = ||(d-^f(a;.Mx)||^ + ||pf(a;„)^(x)||^ 

= \m(x)\\l 1 + \\ P H^)m\\ 2 gM - (4-4) 



The second term of the RHS of Q4.4| ) proves to be expressed as 

\\p x (uj x )ifj(x)\\ 2 gM = g M 1 {p x (uj x )ip(x),p x (uj x )ip(x)) 

= -9m(iI>(z),P%(w x ) ® p$(v x )ij)(x)) 

= (${x),A*${x)), (4.5) 

where we have to make remarks on the notations used; the g^i > ) denotes the inner product on 
T*M © H x , the product p$(u x ) © p x {oj x ) is to be understood as a tensor product in T*M © T*M © 
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EndTt x and the A x is denned by 

A* := -g£{x) o (pX(u x ) ® p*K)) G EndW* (4.6) 

with g^j(x) taken as the inner product on T*M T*M. We notice also that at the second line of 
we have used the fact that (p$(£)v,v') = — (v, p x (£)v') for any £ G g and for any v,v' G Tt x . 



The equivariance of the connection form (3.3) and the invariance of the metric are put together to 
imply that 

A* x = pX( g )AXpX( 9 - 1 ). (4.7) 

Then we observe that the operator A x acting on L2(M;TC X ) through (A x tp)(x) := A x ip(x) leaves 
L2{M]TL X ) G invariant, so that one obtains A x G End(?i x ) We can put the A x in another form. 
Since the inner product gM(%) '■ T X M®T X M — > R gives rise to an isomorphism gM(x) ■ T X M — > T*M, 
its inverse g M l (x) : T^M — > T X M induces an inner product on the cotangent space g~^(x) : T*M <8> 
T*M — > R in the dual manner. Hence g~^ is viewed as a symmetric tensor field gT^ : M — > TM®TM. 
In terms of gj^ 1 along with the connection form cj^ : T X M — > g/gj. and the representation of Lie 
algebra p* : g — > EndW-^, the tensor field A x : M — > EndW x takes the form 

yix :=-(^0^) o ( w0w ) O9 -i. ( 48 ) 

If p x is not a trivial representation and if G acts on M nontrivially, then A x is a positive definite 
operator. We call the A x the rotational energy operator, the reason for which comes from molecular 
mechanics with G = SO(3). In fact, when applied to molecular mechanics, the quantity (4.5) is 
interpreted as the rotational energy density. 

We proceed to study the first term of the RHS of ( |4.4| ). To this end, we have to make a review 
of Hodge's star operator *m '■ f2 k (M) — ► Q m ~ k (M), which is defined through 

a A * M /3 = (a, (3) 9M v M , a, (3 G f2 k (M), (4.9) 

where (a, (3) in the RHS is the inner product on f\ k T*M defined by the metric gu- Let {e%, • • • , e m } 
be a local orthonormal frame field on M with respect to the metric gu- Then the star operator *m 
is explicitly given by 

-i m 

*mo> = — ( i (ej 1 )---i{e jk )a)(i(e jl )---i(e jk )v M )- (4.10) 

' ji,-,jfc=l 

It is easily verified that *m *m ct = (— l) k ( m ~ k ^a. The defining equation ( |4.1| ) of the Laplacian Am 
is then rewritten as 

\\df(x)\\ 2 gM v M = I [ dfA* M df 



M JM 



d(f *m df)- fd * M df 

M JM 
-1. 



f(* M d*M df)v M , (4.11) 

M 

where we have used Stokes' theorem to eliminate the first term on the second line. Then the Laplacian 
takes the form 

Am/ = - *m d *M df = - * M d * M df. (4.12) 
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In terms of local coordinates (x , • • • , x m ) of M, the metric and the volume form are expressed as 
9m = J2i,j 9Mij{x) dx 1 % dxi and vm = vm{%) dx 1 A • • • A dx m = det gMij(x) dx 1 A • • • A dx m , 
respectively. The Laplacian Am is then expressed as 

A«/ = -^E iS l(»M(,)(W)«W^). (4-13) 



as is well-known. 



We examine (4.12) more closely, using the horizontal-vertical decomposition T X M = H x © V x 



The measure \im of M projects to a measure /ig of Q through it : M — > Q; 



<p(q)dfiQ(q) := (<p o n){x)dn M {x), tp G G U (Q). (4.14) 

In what follows, we restrict ourselves to the maximum stratum JkL, which is an open and dense 
subset of M. Let t>g be a volume form on associated with the measure \iq. We can define a 
Riemannian metric <?q on through 7t*<?q = g^lH, where gu\H denote the restriction of the metric 
9m '■ T X M x T X M — > R to the horizontal subspace; 9m\H : H x x iT x — > fi. Note that the definition 
of qq is independent of the choice of x G 7r _1 (7r(x)), because of the G-invariance of gfjvf- The map 
7r : (Mui 9m) (Qli,9q) then becomes a Riemannian submersion. It is to be noted that the volume 
form vq does not coincides with the volume form induced from the metric qq. We denote the set of 
all the horizontal and the vertical fc-forms on by 

Q k H {M^) := {a G Q k {M^) \ i(X)a = 0, VX G V x , Vx G M^}, (4.15) 
ft$(Af„) := {a G fi k {M^) \ i(X)a = 0, VX G ,ff x , G M^}, (4.16) 

respectively. Note that we have already put dimM = m, dimQ M = n = m — p. We define the 
horizontal and the vertical volume forms, vh £ ^h(-^i) ana - V V ^ ^yC^u)j through vh = tt*vq 
and vm = vh h W, respectively. It appears that the forms vh and vy are uniquely determined and 
G- invariant. 

Moreover, it can be shown that both vh and vy are closed forms. It is easy to verify that vh 
is closed; dvn = d(n*VQ) = TT*(dvQ) = 0, since vq is a top form of Q^. To prove that dvy = 0, we 
use a local trivialization over an open set W C Q^; n (W) C is identified with W x F, where 
F := G/Gn is the maximum orbit. The trivialization induces a surjective map np : ir~ 1 (W) — ► -F 
which is G-equivariant, that is, np(gx) = girp(x) for each g £ G and x G 7r _1 (l^ 7 ). A restriction of 
the map 7Tp to each fiber gives a diffeomorphism ir~ 1 (q) = F for each q G W. Let vj? be a G-invariant 
volume form on F which is normalized as J f vf = 1. Then vjr is uniquely determined. Now the 



definition (4.14) of the measure \xq is put in the form 



( P{q)vq = / ((po7r)(x)v M = / {fOTr)(x)v H A vy 

<P(q)vQ[[ vy). (4.17) 

Since tp is arbitrary, this equation implies that the volume form vy restricted to each fiber 7r -1 (g) is 
also normalized as /„■- vy = 1 for each q G Q^, so that vy is a G-invariant normalized volume form 
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on each fiber ir^ 1 (q). It then follows that vy = it f vf- As a consequence, one has dvy = n F (dvF) = 0, 
since vp is a top form of F. 

Using vh and vy, we define the horizontal and the vertical star operators *h '■ ^h(^^) ~~ * 
0]f k (M tl ) and *y : Qy{M^) -> Q 1 v ~ k {M lx ) through 

*M a = (* H a) Avy, a G (A; < n), (4.18) 

*M/S = (-l) n ^A*y/?, /J G ^(M M ) (fc < p), (4.19) 

respectively. According to the decomposition T*M = H*®V*, we break up df into df = {df)jj+{df)y. 
Then *A/d/ is accordingly expressed as 

*M df = * H (df) H A vy + (-l) n ^ A *y(#>- (4.20) 

Since cfo# = and <itv = 0, we obtain 

d* M df= (d * H (df) H ) A vy + v H A (d *y (df)y). (4.21) 

Thus Eq. (|4.12j ) is expressed as 

- A M f = *lld * M df = * H \d * H (df) H ) H + *y 1 (d *y (df)y)y, (4.22) 

which means that the Laplacian Am is decomposed into horizontal and vertical components. 

The above argument can be extended to H x ® Q k (M^) 9* Q k {M^H x ) and to Q k H {M^H x ) G 
straightforwardly; the star operators *m and *h are extended to be applicable to H x -valued forms 
on and to 7^ x -valued horizontal forms on M„, respectively. Hence, for an equivariant function 
ip G C™{M^H X ) G , Eq.flP^) S ives rise to 

- A x ip = * M l d* M dtp = * H 1 (d* H (dlp) H )H + * v 1 {d*y (dlp)y)y. (4.23) 

For tp G C£°(M IM ; Tt x ) G , we have (dip)u = Dip by the definition of the covariant derivation ( p.ll ). 
In view of the first term of the RHS of Q4.23 ), we are led to the definition of the adjoint operator 
Dt : Q k + l {M^H x ) G - Q k H {M^H x ) G of D; 

£>t := - *^ = -(-l) fe (™- fe ) * H D* H . (4.24) 

By using (|4.2|) , ( ^4.4[) , ([4.5|), (|4.22 ), and ( 4.24 ), we accomplish the decomposition of the Laplacian A x 
into the horizontal and the vertical components; 

A x iP = D^DiP + A x ip, ip G C™XM^H X ) G . (4.25) 

We call the D^D the vibrational energy operator for the reason that the integral (4.2) is interpreted 
as the vibrational energy when our general formalism is applied to molecular mechanics. 

As was noted above, the smoothness of the horizontal-vertical decomposition is ensured only in 
the open dense subset C M, so that the RHS of ( [4.25 ) makes sense only in M^. However, A x 
is actually a self-adjoint operator on L2(M;Tt x ) by definition, so that the LHS of (4.2) holds for ip 
defined throughout M. Hence we may expect that some boundary condition arises on dM^ = M s i ng 
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to make D*D + A x into a self-adjoint operator on L2(M;H. X ). Since the requirement imposed on ip 
is that ij) is to be equivariant, we obtain the boundary condition on dM^, 



p x (g)^j(x) = %b(x), g eG x ,xe dM^ (4.26) 

or, in the form of derivative, 

pUZMx) =0, e e g,, x G dMfj,. (4.27) 

In case of G x ^ {e} for the maximum orbit type /i, the equivariant functions are, of course, subject 
to the condition p x (g)i/;(x) = tp(x) for g G G x , x G M^. The condition ( 4.26 ) or ( 4.27] ) states that 



at singular points the equivariant functions are subject to a stronger condition since dimG x rises up 
at singular points. 

C. Angular momentum and inertia tensor 

Now we wish to introduce the angular momentum and the inertia tensor, which are closely related 
with the connection form and the rotational energy operator A x . 

The angular momentum is a map L : T*M — ► g*; which is defined through 

(L x (p),0 := (p,9 x (0), P G T*M, $ G g, (4.28) 

where (-,-)'s in the LHS and in the RHS denote the pairing between g* and g and that between 
T*M and T X M, respectively, and X {£){= £m(x)) is the infinitesimal generator induced by £ G g. 
The angular momentum L : T*M — > g* is a typical example of momentum maps due to Marsden 
and Souriau[26|. By the use of the isomorphism gu ■ TM — > T*M, the angular momentum can be 



rewritten as the map L := L o g~M '■ TM — > g*; u i— ► L x (v), which is expressed as 

<L x (u), := <7Af («, X (Q), v G TjjM, £ G g. (4.29) 

Namely, L is a g*-valued one-form on M, which we call the angular momentum form. 
The inertia tensor is a tensor field I : M — > g* (g> g*; x I x , which is defined through 

4(e, o := 5a/(^(o, uo), e, c e g. (4.30) 

On account of # fl;E (Ad 9 <!;) = g*(0 x {£)) and of 5* #m = ffM for any g £ G, the inertia tensor transforms 
according to 

7 g:B (Ad^,Ad fl C) = 4(C ! C)- (4-31) 

In other word, the map I : M — > g* ® g* is equivariant; 5*/ = (Ad*_i (g) Ad*_i)J. For an arbitrary 
x & M fixed, the quadratic form I x : g <S> g — > -R can be regarded as a map I x : g — ► g* ; which 
has Ker I x = g x and Iml x = {0 G g* | = 0, V£ G g^,} = (g/g x )*. Then it can give rise to an 

isomorphism I x : g/g x —*(g/s x )*' Hence, there exists the inverse {I x )~ l ■ (g/Sx)* ^g/Sx; which is 
identified with a quadratic form (I x ) ■ (g/g x .)* ® (s/SxT ~~ ^ The map / x will be referred to as 
an inertia operator. The inertia operator is called the locked inertia tensor by Simo et al.[^7]]. The 
inertia operator was first introduced by Guichardet|J, and used in || f?J to break up the total energy 
into the sum of rotational and vibrational energies. 
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Prom definition, the angular momentum, the inertia tensor, and the connection form turn out to 
be related by 

L = I o cj, 



where each symbol is to be understood as follows: L x : T X M — > g*, tu x 
Ix '■ e/e x ~^ E*- A proof of (|4.32|) runs as follows: First, the identity 



T X M 



l o e = i, 

can be proved by a straightforward calculation. In fact, for any £, £ G g, we have 

((Lo8 x )(0X)=gM(e x (0,e x (()) = Ut,() = (U0X). 

Next, from the identity ( |4.33| ) it is deduced that 

Lo8ouj = Iouj. 



(4.32) 

(4.33) 
(4.34) 
(4.35) 



From the definition of the angular momentum form ( 4.2S| ), it can be shown that Ker L x = H x , so 
that Eq.( 4.32| ) holds on H x . Moreover, since 9 X o uj x : T X M = V x © H x — > V x is a projection as 
was noted at (p 7 ^), Eq.( |4.35| ) shows that ( 4.32| ) holds on V x . The proof is thus completed. Since 
lmL x = (g/g x )*, we may rewrite ( 4.32| ) as L = I o uj, so that the connection form is put in the form 



cj = I- L oL. 



(4.36) 



This formula will be used to write out the connection form u> in molecular mechanics. 

Owing to the definition of the inertia tensor, I x = gu (Ox <8> Ox), the metric givi\V restricted to 
the vertical subspace takes the form 



9m\V = I o (u> (8) u) = (L, cj) , 



(4.37) 



where use has been made of X o lo x \V x = idy x , another expression of (p\4|). Equation ( 4.37 ) can be 
looked upon as describing the rotational energy in classical mechanics. Since the reduced quadratic 
form I x = gu ° {0 X <8> X ) : e/s x ® e/e x — > -R is nondegenerate, it has the inverse quadratic form 
(/a;) -1 = g~^l o (u* ®tu x ) : (e/Sx)* ® {e/ExT ~^ which is expressed as a tensor field, x i— ► (7 a ) _1 G 

{I x )- X = (lu x oj x ) o gll (4.38) 

where 5^ is regarded as a symmetric tensor field M — > TM ®TM . From ( |4.8| ) and (4.38), we obtain 
a formula to express A x , in terms of the inertia tensor, as 



A* = -(p*®pX)or\ 



(4.39) 



D. Coordinate representation 



Now we wish to provide a coordinate representation of the reduced Laplacian given in ( 4.25 ). Take 
local coordinates (g 1 , q 2 , • • • , q n ) on an open subset W of Q^. Then the metric gQ and the volume 
form vq take the form gQ = YUj 9Qij(o) ^ ® dq 3 and vq = VQ(q) dq 1 A • • • A dq n , respectively. For 
each x £ tt^ 1 (W) C M^, we take a basis {£i(x), • • • ,£, p (x)} of g/g a , where we have used the same 
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notation for elements of g/g x as those for g for simplicity. Then the components of the reduced 
inertia tensor I x are defined by 

(Tx)af}'-=Tx{Sa(x),£ (x)), a, f3 = 1, • ■ • ,p = dimg/ g;E , (4.40) 

which gives a symmetric positive definite matrix of rank p. The components of its inverse are denoted 
by (I~ 1 ) al3 . Let a be a local section on W of the bundle 7r : — ► Q^, that is, a differentiable map 
a : W — » M M such that 7r o a = idw- Then an equivariant function ■0 G C£°(M; H X ) G is pulled back 
to a function cr*^ : — > W x of the coordinates (g 1 , <? 2 , ■ ■ ■ ,q n ). Our aim is to obtain a coordinate 
expression of a*A x ip, according to the decomposition A x = Z^D + A x . The covariant derivative 
( 3. 13| ) with a = ip is expressed, in terms of (g , g 2 , • • • , g n ), as 

(<r*Ity)( 9 ) = E(^-(P^° (^)) {°*1>){*W- ( 4 - 41 ) 



It then turns out from Q4.24Q , ( j4.25j ), and Q4.39Q that 



VQmjnq) (Jj - Wf ° " o a*) (^-)) 
" E (5"(l))^Wf ° o Wf o o ff )( g )(^)(?). (4.42) 

a,/3=l 

It should be noted here that uj(a*(d/dq 1 )) and p x (uj{a 1f {d / dq 1 ))) denote a component of a "gauge 
potential" and its representation as a matrix acting on 7i x , respectively. Further, p x (£ a (cr(q))) denote 
a matrix representation of the infinitesimal generator (or the angular momentum operator) induced 
by (, a (x) S g. This is one of our main results. 

If there exists a global section a : — ► M„, the fiber bundle ir : M« — > becomes a trivial 
bundle; M M = Q M x (G/G^), and cr*V' becomes a smooth function over the entire domain Q^. In this 
case, the procedure of reduction is nothing but separation of variables. The reduction method is a 
generalization of the method of separation of variables. 

Equation ( |4.42|) makes sense only in the maximum component, Q^, of the orbit space Q. At a 
singular point q G dQ^, the rank of the inertia tensor I a {q) decreases abruptly. As was noticed earlier, 
the equivariance condition provides the boundary condition imposed on ijj and hence on a*tp, which 
is put in the form 

pHO(^)(q) = o, Ceg^). (4.43) 

If Ea(q) = {0} m this imposes no condition on a*ip. At singular points q 6 dQ^, the dimension 
of the isotropy algebra g a r q \ jumps up, so that the value of a*Tp is more strongly restricted there. 
The operator ( [4.42j ) should be accompanied by the condition ( 4.43 ). This is another one of our main 
results. 

E. Example: quantum mechanics on a plane 

Quantum mechanics in a two-dimensional Euclidean space provides a simple but nontrivial example 
of the formulation constructed above. 
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First we take M = R 2 and G = SO (2). Then 50(2) acts on R 2 in the usual manner. Let (r, (f>) 
be the polar coordinate of R 2 . Then the orbit space becomes Q = R>o, which is nothing but the 
radius coordinate r > 0. A point with has a maximum orbit S . The origin r = has a singular 



orbit, that consists of a single point {0}. In this example we have M u 



R 2 



{0}, dM u = {0}, 



Q u = R >0 , and dQ u = {0}. Since R 2 - {0} = R >0 x 5 1 , the fiber bundle it : — > Q u is trivial. 
Then the reduction method becomes the method of separation of variables in this case. 

The space R 2 is equipped with the standard metric gu = dr 2 + r 2 d(j) 2 and the standard measure 
dp-M = rdrdcj). Then the projected measure of R>o is dfiQ = 2nrdr. At each maximum point x G R 2 



with r / 0, the vertical and horizontal subspaces are given by 



V, 



Hi, 



Hi 



R- 

dr ' 



(4.44) 



R . The metric on .R>o is 
— > Qp, . It should be noted 



respectively. At the singular point r = 0, one has Vb = {0} and i^o - 
5q = dr 2 , and thereby one obtains the Riemannian submersion it : M t 
that dfiQ = 2nrdr does not coincide with the metric volume dr. 

Any irreducible unitary representation of 5*0(2) is one-dimensional. It is labeled by an integer 
n G Z and defined by 



Pn : SO{2) -> i7(l); 



cos a 
sin a 



- sin a 
cos a 



(4.45) 



Accordingly, an equivariant function ip : R 2 — > C which satisfies tp(gx) = p n (g)ip(x) becomes a 
function ip n (r,(/>) subject to the condition tp n (r, <j) + a) = e ma ^ n (r, 4>). Thus we can put tp n in the 
form ip n (r, <fi) = e m< ^/ n (r). Here the decomposition ( |2.10[ ) with the projection operators (2.7) realizes 
the ordinary Fourier expansion in angular coordinate, ip(r, 9) = Ylm=-oo ein ^fn{ r )- Since 50(2) is an 
isotropy group at the origin r = 0, smooth equivariant functions ip n must satisfy 



9 I 



0. 



r=0 



(4.46) 
(4.47) 



which has an alternative expression 

nf n {0) = 0. 

This boundary condition illustrates the general condition ( f4,43| ). 

We proceed to reduce the ordinary Laplacian. The metric on the cotangent bundle T*R 2 is 
expressed as 



9m 











d Id 

dr dr r 2 d(j) d(p 



(4.48) 



To obtain the reduced Laplacian, we calculate the integral ( p~2j) for <ip v 



R 2 



\\d^n\\ 2 g M dp M 




1 



rdrd(j) 



+ ^V|/„| 2 }2TTrdr 



j- _l± r dh ^ 

.In \ j ' j t 9 Jn 
r dr dr r z 



2irrdr. 



(4.49) 
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df 

Since f n {r) is bounded as r —* from (14.471), and since r— — — * as r — * from the second line 

dr 



of ( 4.4S| ), the boundary term at r = vanishes. The other term at r = oo vanishes because of the 



assumption that ip n has a compact support. Thus we are left with the reduced Laplacian 



*nfn = ~~ + ^fn- ( 4 - 50 ) 

r dr dr r z 



We note here that the boundary condition (|4.47|) says that / n (0) = for n / and that /o(0) is 



bounded for n = 0. This result gives an example of the general formula ( 4.42Q . An eigenfunction 



of the Laplacian A n with an eigenvalue E > is the ra-th Bessel function J n {y/Er). The Neumann 



function Y n (y/Er) does not satisfy the boundary condition ( 4.47 ). 
F. Rigid body 

We have to notice here that our theory covers quantum mechanics of rigid bodies. In mechanics, 
a rigid body is defined as a collection of mass points in three dimensions, in which their mutual 
distances are kept fixed. In this case, the symmetry group G is 50(3) and the configuration space 
M becomes a single G-orbit. Hence the orbit space Q reduces to a single point. When the orbit is 
maximal, M becomes SO (3) and the inertia tensor I is nondegenerated. When the orbit is singular, 
M is isomorphic to S 2 or a single point, and the inertia tensor is of rank two or zero, respectively. For 
a rigid body, the horizontal component of the Laplacian ( 4.25| ), or the vibrational energy operator, 



vanishes, and therefore the Laplacian reduces to the rotational energy operator; 

A x = A x = -(p x ®p x )or 1 , (4.51) 

which is the Casimir operator acting on the representation space 7i x up to a normalization constant. 
In the language of physics, since the rigid body executes no vibrational motion, it has only rotational 
energy, which is determined by the angular momentum. A simple example will be given in Sec. V. 

V. Quantum molecular mechanics 
A. Jacobi vectors 

In the previous sections, we have set up a general formulation for reduction of quantum dynamical 
systems on the configuration space M with symmetry G. The Hilbert space L2(M) is decomposed 
into the orthogonal direct sum of the spaces of equivariant functions according to the irreducible 



unitary representations of G, as was shown in (2.35). The Laplacian Am is accordingly reduced 



to the operator A x of fl4.25| ) acting on each space of equivariant functions L2(M;7i x ) G . We have 
studied quantum mechanics on M = R 2 with symmetry G = SO (2) to give a concrete example. It 
was a well-known but nontrivial example in which reduction by symmetry serves as the method of 
separating variables. 

Here, we wish to apply the general formulation to molecular mechanics, which is the original 
problem that motivates us. We consider a molecular system consisting of iV atoms in R?. The 
configuration of the molecule is described as an ennuple (xi, ■ ■ ■ , jbjv) °f the position, ajj G R?, of 
each atom. Masses of the atoms are denoted by (mi, • • • , mjv) with mi G R>o- Assume that we are 
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working with the center-of-mass system; 

N 

M = {(xi, ■ ■ ■ , x N ) G (R 3 ) N | m i x i = 0}, (5-1) 

i=i 

which is a linear subspace of (R 3 ) N . Let g G G = 50(3) act on x = (aci, • • • , xjv) £ M by 
= (paii, • • • ,gxw). We call M and Q = M/SO(3) the molecular configuration space and the shape 
space, respectively. We may regard x = (x\, ■ ■ ■ ,xn) as a 3 x N matrix. According as the rank 
of x is 3, 2, 1, or 0, the configuration x is called a generic configuration, a planar one, a collinear 
one, or a collision one, respectively. A generic or planar configuration has a maximum orbit which 
is diffeomorphic to SO (3). A collinear configuration, in which all the atoms are placed along a line, 
has a singular orbit which is diffeomorphic to S 2 . A collision configuration x = (0, • • • , 0) has another 
singular orbit which consists of a single point. We are going to review a geometric setting on M in 
what follows. The topology of the shape space Q will be studied in the next subsection. 

While a tangent vector v G T X M is denoted by v = (v\, • • • , vn) G (R 3 ) N along with the condition 
J2iLi mi v i = 0, a cotangent vector p G T*M is denoted by p = (p 1 , ■ ■ ■ ,Pjy) G (R 3 ) N along with the 
condition J2iL\Pi = 0. The pairing between T X M and T*M is given by (p, v) := X)i=i(Pj> v i)i where 
(•, •) denotes the standard inner product of i? 3 . Each element of the Lie algebra £ G so(3) is identified 
with a vector £ G R 3 and induces the infinitesimal transformation 9 X (£) = (£ x x\, ■■■ ,£ x ajjv), 
where x means the standard vector product in R s . Moreover, the dual space of the Lie algebra 
so (3)* is also identified with R 3 . 

The space (R' i ) N is equipped with a Riemannian metric K; 

N 

K :=^2 m i(dxi, dxi). (5.2) 
i=l 

Although K is twice the kinetic energy, we will call it the kinetic energy simply. The subspace 
M C (R?) N inherits the Riemannian metric K. From the definition ( [4.28 ), the angular momentum 
L : T*M -> so(3)* = R 3 is expressed as 

N N 

{L x (p),i) = (p,9 x (£)}, = J2(Pvt x x i) = (E x * x Pi'O, (5-3) 

1=1 8=1 

and hence takes the usual form 

N 



L = Y,XiX Pi . (5.4) 



=i 



According to (4.29), the angular momentum form L : TM — » so(3)* is expressed as 

N N 

(L x (v),£) = K{v,0 x {£)) = ^2mi{vi,£, x x { ) = (^rriiXi x (5.5) 

i=l i=l 

which implies that 

iV 

L = mjCEj x dajj. (5-6) 

i=l 

The Jacobi vectors are of great help in describing many-body systems. Let us remind us of the 
definition of the Jacobi vectors. By Mi G -R>o and X{ G .R 3 we denote the sum of the masses from 
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the first to the i-th atom and the center of mass of the set of i atoms, respectively, 

i 1 * 

Mi-^mj, X i :=—Y,m j x j (i = 1, • • • , TV). (5.7) 

In particular, one has X\ = xi, and Xn is equal to the center of mass of the whole system. Then 
the Jacobi vectors (r^\r[ N \ • • • , rj^2 1 ) are defined by 

Of course, in the center-of-mass system, one has Xjy = = 0. The Jacobi vectors (r[ N \ ■ ■ • rj^2 1 ) 
provide a coordinate system to M, and give rise to the isomorphism M = (.R 3 )^ -1 ). 

Good use is made of the Jacobi vectors to prove the additivity of the kinetic energy K and of the 
angular momentum L in the number of atoms. In fact, one can verify that 

N N-l 

da*) = £ (drf ),dr{^), (5.9) 

i=l i=0 

L (N) := x dXi = rf ] x drfK (5.10) 

i=l i=0 

The additivity of can be proved by induction with respect to N. A straightforward calculation 
yields 

K (N+i) _ K (N) 

= II^IP + II^II 2 -!!^!! 2 



|2 



= M N+1 \\dX N+1 \\ 2 + I | da . _ d x w || 2 - M^HdXjvl 

Mat + mAr+i 

= —!—\\M N dX N + m N+1 dx N+1 \\ 2 + M ^ mjV+1 ||da;jv+i - dX N \\ 2 - M N \\dX N \\ 2 
+ MjvmAr+i - MjyMAf+1 2 m^ +1 + Mjvmjv+i 2 

= iv^v+i l|d JV " + — %; l|cte7V+l11 

= mAT+iUdajAr+iH 2 . (5.11) 
In a similar manner, the additivity of the angular momentum is verified as follows: 

= +1 > x dr <" +1 > + x cir^ - r ™ x dr™ 

= Mn+iX n+i x dX^ + i 

+ „ , NUlN+1 (x n+1 - X N ) x (dx N+1 - dX N ) - M N X N x dX N 
M N + mjv+i 

(M N X N + mjv+iaJjv+i) x (M N dX N + mjv+icteAr+i) 



+ M f, mjV+1 (a;iv+i - X N ) x (dajiv+i - d^jv) - M N X N x dXjv 

Mjv+1 

M% + M N m N+1 - M N M N+1 m 2 N+l + M N m N+1 

— 2 — Xat x dX^v + ^— T7 jbjv+i x dx N+1 

M N+ i M N+ i 

m N+1 x N+1 x dx N+1 . (5.12) 
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In the following, we fix the number of atoms iV and suppress the superscript (N). 

According to the relations ( }4.33| ) and ( 5.1C| ), the inertia operator I x : so(3) — > so(3)* is defined 
for each £ G so(3) = R 3 by 

JV-l N-l 

= L(0 X (£))= Y / r l x(Hxr l )= £ ((r^K - forjri). (5.13) 

i=0 i=0 

According as x is a generic configuration, a planar one, a collinear one, or the collision one, the rank 
of I x is 3, 3, 2, or 0, respectively. Unfortunately, there is no concise expression for the inverse (I x ) 
in general. However, the connection form lo is expressed, from (4.36), as 



JV-l 

u> = r 1 J2 r i x dr i ( 5 - 14 ) 

To formulate molecular mechanics, we need the invariant volume form vm of M associated to the 
metric K, 

v M = d 3 ri A • • • A d 3 r N ^i. (5.15) 

Thus we have made a geometric setting to apply our formalism to molecular mechanics. However, 
before application we have to examine the topology of the orbit space (or shape space) Q = M/G 
for N- atomic molecules. 

B. Topology of the shape space 

Let M(rn, n) and M(m, n)t- be the vector space ofraxn matrices over R, and the set ofraxn matrices 
of rank k, respectively. By S(n) we denote the set of all the positive semi-definite symmetric n x n 
matrices, and set S(n)i- := S(n) n M(n, n)&. Of course, k < m,n. Let 0(n) denote the orthogonal 
group acting on R n as usual. Then O(m) and 0(n) act on M(m,n) to the left and to the right, 
respectively. We can verify now that 

M ( ™, „), 2 Q(m)/0(m-IQxSW t xO(n-*)\0(,,) ^ 

from the observation of the fact that an arbitrary linear map (p : R n — > i? m of rank can be 
expressed as a composition 99 = 200-071" of three linear maps, where ir, a, and z are an orthogonal 
submersion it : R n — » a positive-definite symmetric operator cr : — > and an orthogonal 
immersion i : R k — > respectively. Here we call a linear map 7r : i? n — ► an orthogonal 

submersion, when it is surjective and satisfies tt o *7r = id on Similarly, we call i : R k ^ m 
an orthogonal immersion, when it is injective and satisfies *i o i = id on R k . The set of all the 
orthogonal submersions {tt : R n —> R k } is identified with a Stiefel manifold 0(n — k)\0(n), while 
the set of all the orthogonal immersions {i : R k — > i? m } is identified with another Stiefel manifold 
0(m) /0(m — k) . An equivalence relation ~ is defined on the triplet (i, a, tt) by the action of g £ 0(k) 
through (i,a,ir) ~ (ig" 1 , g<rg~ l , gir). In particular, for m = n = k, Eq.( |5.16 ) becomes 



GL(n) = M(n, n) n ^ 5(n)„ x O(n), (5.17) 
which is nothing but the so-called polar decomposition of regular matrices. 
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Thus the configuration space of the molecule, M = (R 3 ) N 1 = M(3,N — 1), is identified with 

mp.jv-i)^ U oP)M»-*)^W*»-o^-'-*)\o^-'), (5 . 18) 

0<fc<min(3,7V-l) ^ ' 

Each component with k = 0, 1,2, 3 corresponds to the set of collision, collinear, planar, and generic 
configurations, respectively. Note that a point x £ M is of the maximum type or of the singular 
type, according as k = 2, 3 or k = 0, 1. Strata of the shape space Q = M/SO(3) are then given by 



lW „ 50(3)\0(3)/0(3 - *) x 5(fc) ifc x 0(iV - 1 - fc)\0(JV - 1) 

O(k) 



The topology of strata for few-body problems with N = 3 and TV = 4 is already studied by one of 
the authors^ and Narasimhan-RamadaspSl], respectively. Coordinates of the iV-body problem are 



also studied by Littlejohn and Reinsch[29|. They also wrote a comprehensive review [pQ] on gauge 
fields in the iV-body problem, and studied also complexes of rigid molecules [[31]]. 

We write out topology with TV = 2, 3, 4 to give definite examples. We denote by i?>o the positive 
real numbers (0, oo) and by D n , S n , and RP n an n-dimensional disk, sphere, projective space, 
respectively. 



N 



N = 3: 



* 5(1)! * R >0 , 
Qo ] = {0}. (5.20) 



Q { 2 ] * S(2) z = R >0 x D 2 = R >0 x R 2 , 

Q? = s(i)i x om^om ~ 5(2)1 ~ ^ >0 x Rpl ~ R>0 xSl - R2 ~ {°h 

Q { o ] = {0}. (5.21) 



N = 4: 

Qf = x 5(3) 3 = ^2 x i* >0 x D 5 = R >0 x (5 5 - S 4 ), 

Q? * S(2h X ( 2 ( ; A ° (3) ^ *(3), ^ H>o x - RP% 

g[, 4) = {0}. (5.22) 

In the case of ./V = 3, the union of , fe = 0, 1,2, forms the shape space Q = R 2 x -R>o- The 
maximum stratum is Q% — R 2 x i?>o- For N = 4, the union U ^ s ^ ne maximum stratum, 
which is diffeomorphic to i?>o x (S" 5 — RP 2 ). 
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C. Triatomic molecules 



To make a practical application of the above general formalism, we concentrate on the triatomic 
molecules in the rest of the paper. The configuration space then becomes M = (i? 3 ) 2 = {(ri,r 2 )}. 
The maximum stratum is diffeomorphic with M(3, 2) 2 , the space of 3 x 2 matrices of maximal 
rank. The stratum of singular orbit type, c?M„, is the union M(3, 2)i U M(3, 2)o- Dragt[32| and his 
successors have introduced a useful coordinate system (a, (3, 7, p, x, <t>) of M by setting 

X 4> X 4> 

n = p (cos - cos -u 3 + sin - sin ^u 2 ), (5.23) 

r 2 = p (cos - sin -1*3 - sin - cos - u 2 ). (5.24) 
Here (1*1, 1*2,1*3) is an orthonormal basis of i? 3 parametrized by the Euler angles (a, (3, 7) as 

(tii,U2,ti3) = 0(ei,e 2 ,c 3 ), g = e aj3 e? Ja eY J *, (5.25) 

where (ei, 62,63) is the standard basis of -R 3 and g is an element of SO(3) with (Ji, J2, J3) the 
standard basis of so(3) defined by J{V = ej X v (i = 1,2,3) for each v £ i? 3 . We notice here that 
the orientation of the frame (1*1,1*2,1*3) is different from that of the original article. We choose 
the orientation to bring the collinear configurations into the direction of 1*3. The ranges of the 
coordinates are given by 

< a < 2tt, < (3 < vr, < 7 < 2vr, 

0<p, 0<x<~ 0<^<2vr. (5.26) 

The geometric meaning of (/?, x> <j>) is clear on introducing coordinates (gi, (72, 93) by 

gi := ||ri|| 2 — 1 1 -r*2 1 1 2 = p 2 cosxcos</>, (5-27) 
52 := 2(r-!,r 2 ) = p 2 cos x sin 0, (5.28) 
g 3 := 2||n x r 2 || = p 2 sin^. (5.29) 

They are invariant under the action of S*0(3) on M, and provide the projection ir : M — > Q = 
M/SO(3); (n, r 2 ) 1— > (gi, 52, 93)- With this expression of 7T, it is easy to show that Q is homeomorphic 
to the upper half space _R| = i? 2 x i?>Q. The space, 7r(M M ), of maximum orbits is diffeomorphic with 
R 2 x R>o- The boundary surface, R 2 x {0}, determined by g 3 = or x = describes the orbit space 
for collinear configurations, in which the atoms make a line along 1*3. The origin (0,0,0) represents 
the collision configuration. These observations coincide with ( 5.21 ). On the other hand, the Euler 
angles (a, (3, 7) are regarded as a coordinate system of the fiber of the bundle ir : M M — > Q^. Note 
that, in the set of singular points dM^, one has T\ = pcos — 1*3, r 2 = psin — 1*3 with 1*3 = e aj3 e^ j2 e 3 , 
which shows that (a,/?) serve as coordinates for the orbit diffeomorphic with S 2 , and that (p,<j>) are 
coordinates for dQ^ = R 2 . 

In terms of the coordinates introduced above, we are to write out the geometric objects in the 
explicit form. It is also convenient for later use to introduce the Maurer-Cartan one-form; 

9~ X dg 

= J\ (sin 7 d(3 — sin [3 cos 7 da) + J 2 (cos 7 d/3 + sin (3 sin 7 da) + J 3 ( d*y + cos /3 da) 
= Ji©i + J 2 <9 2 + J 3 6> 3 - (5.30) 
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Then the metric g\j = K obtained in (5.9) takes the form 
9M = ||dri|| 2 + ||dr 2 || 2 

i 2 ^ 2 t 2 1 2 7/2 

= dp + -p dx + -p cos x dp 



+p 2 (©i - - sin x 



2 , 2 2 X ^,2 , 2 ■ 2 X ^,2 

+ p cos 2&2 + P sm 2 3- 



The angular momentum form L : TM — ► so(3)* = i? 3 obtained in ( |5.10| ) is expressed as 



n x dn + r 2 x 
1 

2 



p 2 wi (@i sin x dtp) + p 2 cos 2 — ui 02 + p 2 sin 2 — U3 ©3. 



2 ' 2 

The inertia operator I x : so(3) = ii 3 — ► so(3)* = R 3 obtained in ( 5.13j ) is put in the form 
4(£) = ri x (£ x n) + r 2 x (£ x r 2 ) 

2 



p 2 (^,Mi)ui + p 2 cos 2 %(£,u 2 ) u 2 + p 2 sin 2 -(£,1*3) **3- 



The connection form cj : TM — > so(3) = i? is then written out, according to ( 4.36 ), as 

w = I~ l L = t*i ((9i - - sinx#) + u 2 02 + u 3 3 . 
Hence the vertical component of the metric ( 4.37| ) is given by 



9m\V = I o (it) <g> u>) = (L, u) 



P 2 (°i ~ 2 s[nx 



1 2 2 A ^2 1 2 ■ 2A n 
+ p COS - @ 2 + p sm - 3 



and thereby the metric gQ such that ir*gQ = Qm\H becomes 

9Q = dp 2 + ^p 2 dx 2 + ^p 2 cos 2 x < 

Then its inverse is a tensor field given by 

, , , d d 4 d d A 
(sq) = ir®JT + -2 IT ® IT + 



9p 5p p 2 Ox dx p 2 cos 2 x 
The volume form defined by ( |5.31 ) is of the form 



d d 

777 ® 7T7 



1 

16 



p 5 sin 2x dp A dx A d(f> A 0i A 2 A 03- 



Then the volume form associated with the measure ( |4.14| ) is 



(5.31) 



(5.32) 



(5.33) 



(5.34) 



(5.35) 



(5.36) 



(5.37) 



(5.38) 



7T" 



V Q = Y P Sin 2 ^ d pAdxA~ 



(5.39) 



since / 6>i A 2 A ©3 = 8vr 2 . 

To describe quantum mechanics for the triatomic molecule, we need the Hilbert space of sections 
of vector bundles associated with the stratification. Any irreducible unitary representation of 5*0(3) 
is characterized by a nonneg ative integer i and denoted by p l : 50(3) -> U{C 2i+1 ). We put J- := 
(//)*(ej) for the standard basis G i? 3 = so(3). Since the base space of the stratified fiber bundle 
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7r : M —> Q = R x R>o is contractible, the bundle is a trivial bundle. Through a global section 
a : Q -> M with a = /3 = 7 = 0in ( |5.23|) , ( |5.24[) , and ( |5.25[) , any equivariant function ip : M C 2m 
is pulled back to a function & := a* ij) : i?> 
form 

J 3 # = for x = 



C 2 ^ +1 . Then the boundary condition ( [4.43 ) takes the 

(5.40) 



on M(3,2)i and 



J x \p = J 2 xp = J xjr = o for 



p = (5.41) 

on M(3,2) , respectively, with the coordinate system defined at (5.27), ( |5.28| ), and ( |5.29| ). What 
Eq. (pilcD means is as follows: At a collinear configuration determined by x = 0, the molecule 
lying in the line along 113 has the vanishing angular momentum about u^; (p )* (1*3)^(2:) = 0, so 
that one has p l {g)(p e )*{es)p e (g~ 1 )ip(ga(q)) = with 1*3 = ge^, which is equivalent to ( 5.40| ). Since 
J3 = diag(^, £ — 1, ... ,0, ... , —£), the components \F m of \P with m/0 vanish, if £ / 0. Furthermore, 
Eq. ( p"4lD means that at the collision configuration determined by p = 0, the molecule cannot 
carry nonzero angular momentum, so that the wave function can have a nonvanishing value only 
when £ = 0. These conditions are analogs of that for the two-dimensional case ( [4.47| ). Now we have 
implemented the consideration of singular case by providing the boundary condition ( ^.40|) and Q5.41 ) 
which we skipped in the previous workQ. 

In conclusion, we write down the Laplacian flP^ ), combining flOaD , ( fjjj ), O), and O); 

-£&{p,xA) 

1 f5 5 . n 9 3 4 5 . 3 

= 5 ■ o 1 sm2x— + -p sin2x— 

p° sin Z\ I "X P ®X 

d 1 o\ 4 



+ - sin xJi 
1 



2 * 7 p 2 cos 2 x 



-~2 i (Ji) 2 H 2? — ToT^ 2 

p 2 cos 2 (x/2) 



+ 



p sin2x 
1 



5 



■ sin x«/i 



y 2 



5 d 4 
<9p 2 p dp p 2 I <9x 2 



1 



+ 2 cot 2x 



sin 2 (x/2 
d\ ! 

<9x, 



+ 



p 2 cos 2 x 



1.- 

9^ + 2 SmxJl 



p^ I cos / (x/2) 



+ 



1 



sin 2 (x/2; 



■(Ja 



(5.42) 



which reproduces the result of 0] . The first and last terms of the RHS of ( |5.42| ) are vibrational and 
rotational energy operators, up to sign, respectively. As was pointed out in [[?]], if the vibrational 
energy operator is separated off, and if the internal coordinates (p, x> <P) are fixed, the operator |A 
reduces to the well known Hamiltonian for a rigid rotor of plane body. 



VI. Concluding remarks 

In this paper we formulated the general method of reduction of quantum systems with symmetry by 
the use of the Peter- Weyl theorem. Although the method is well-known implicitly among Physicists, 
we developed it explicitly to give rigorous grounds to quantum mechanics describing molecular mo- 
tions. We studied the stratification of manifolds according to the action of a symmetry Lie group 
and then defined a stratified bundle and a stratified connection as generalization of fiber bundles and 
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connections. Further, we showed that the reduced quantum system is a pair of the Hilbert space 
and the Hamiltonian which are the space of equivariant functions and the Laplacian expressed in 
terms of covariant derivation with the stratified connection, respectively. We found the boundary 
condition that is imposed on the equivariant functions to make the reduced Laplacian a self-adjoint 
operator. Finally, the general formulation for reduction was applied to iV-atomic molecules, and 
triatomic molecules were examined in particular. 

The stratified connection on the stratified bundle is newly introduced as a generalization of 
connections on principal fiber bundles and is used to describe the reduced Laplacian. One of our main 
results is to have determined the boundary condition for making the reduced Laplacian self-adjoint. 
Emmrich and H. R6mer|33|] analyzed Laplacians on orbifolds to study quantization of systems with 
gauge symmetry. They found that the Laplacian on an orbifold is not essentially self- adjoint and 
therefore its self-adjoint extension is not unique. According to our method, the reduced Laplacian 
is self-adjoint by its definition and the boundary condition is accompanied automatically by the 
symmetry. 

We would like to make some comments on remaining problems. First, although we built a 
general formulation to deal with quantum molecular systems, we do not obtain spectrum of the 
reduced Hamiltonian. Even the three-body problem is difficult to solve. It is desired to develop an 
approximate method to solve the reduced eigenvalue problem of physically interesting systems. 

Second, for application to real molecules, electronic structure must be considered. Of course, 
spins and statistics of electrons and nucleus also must be taken into account. These are left for 
future investigation. 
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